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Abstract 

We study a random graph model which is a superposition of the bond percolation 
model on Z d with probability p of an edge, and a classical random graph G(n,c/n). 
We show that this model, being a homogeneous random graph, has a natural relation 
to the so-called "rank 1 case" of inhomogeneous random graphs. This allows us to use 
the newly developed theory of inhomogeneous random graphs to describe the phase 
diagram on the set of parameters c > and < p < p c , where p c = p c (d) is the 
critical probability for the bond percolation on Z d . The phase transition is similar to 
the classical random graph, it is of the second order. We also find the scaled size of 
the largest connected component above the phase transition. 

1 Introduction. 

We consider a graph on the set of vertices B(N) := {— N, . . . , N} d in Z d , d > 1, with two 
types of edges: the short-range edges connect independently with probability p each pair u 
and v if \u — v\ = 1, and the long-range edges connect independently any pair of two vertices 
with probability c/\B(N)\. (Here for any set A we denote \A\ the number of the elements 
in A.) This graph, call it Gn(p,c) is a superposition of the bond percolation model (see, 
e.g., [4]), where each pair of neighbours in Z d is connected with probability p, and a random 
graph model G nfi / n (see, e.g., [6]) on n vertices, where each vertex is connected to any other 
vertex with probability c/n; all the edges in both models are independent. By this definition 
there can be one or two edges between two vertices in graph Gn(p,c), and in the last case 
the edges are of different types. 

The introduced model is a simplification of the most common graphs designed to study 
natural phenomena, in particular, biological neural networks [12]. Notice the difference 
between Gn(p,c) and the so-called " small- world" models intensively studied after [13]. In 
the " small- world" models where edges from the grid may be kept or removed, only finite 
number (often at most 2d) of the long-range edges may come out of each vertex, and the 
probability of those is a fixed number. 
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We are interested in the connectivity of the introduced graph G^(p, c) as iV — > oo. We 
say that two vertices are connected, if there is a path of edges, no matter of which types, 
between them. Clearly, if c = 0, we have a purely bond percolation model on Z d , where any 
edge from the grid is kept (i.e., "is open" in the terminology of percolation theory) with a 
probability p, or, alternatively, removed with a probability 1 — p. Let us recall some basic 
facts from the percolation theory which we need here. Denote C an open cluster containing 
the origin of Z d in the bond percolation model. It is known (see, e.g., [3]) that for any d > 1 
there is p c = p c (d) such that 



P{|C| = oo} 



= 0, if p < p c 
> 0, if p > p c 



where < p c < 1, unless d — 1, in which case, obviously, p c — 1. We shall assume here that 
< p < p c , in which case the connected components formed by the short-range edges only, 
are finite with probability one. Recall also that for all < p < p c the limit 

C(p) = limn^ \^~\ogP {\C\ = n}j (1.1) 

exists and satisfies ((p) > (Theorem (6.78) from [1]). 

Let further C\ {g^ denote the size (the number of vertices) of the largest connected 
component in a graph G. 

Theorem 1.1. Assume, that d> 1 and < p < p c {d). Define 

c CT ( P ) = -^j. (1.2) 

i) If c < c cr (p) define y to be the root of E c\C\e c ^ c ^ v = 1, and set 

a(p,c) := (c + ay-Ece^y 1 . (1.3) 

Then for any a > a(p, c) 

P{C 1 (G JV (P, C )) >alog| J B(iV)|} ^ 0. (1.4) 

as N — >• oo. 

ii) If c> c cr (p) then 

Wn)\ U p (L5) 



2 



as N — > oo ; with (3 = /3(p, c) defined as the maximal solution to 

(3 = 1-E{e- C/3|C| }. (1.6) 

In view of (11.11) it is obvious that xip) '■= E|C| < oo for all < p < p c . It is also known 
(see Theorem (6.108) and (6.52) in [lj) that x{p) is analytic function of p on [0,p c ) and 
x(p) — > oo as p — > p c . This implies that c cr (p) is continuous, strictly decreasing function on 
[0,p c ) with c cr (0) = 1 and c cr (p c ) = 0. Hence, c cr has inverse, i.e., for any < c < 1 there 
is a unique < p cr (c) < p c = p c {d) such that c cr (p cr (c)) = c. This leads to the following 
duality of the result in Theorem IX . 11 

Corollary 1.1. For any < c < 1 there is a unique < p cr (c) < p c such that for any 
p cr (c) < p < p c graph Gn(p, c) has a giant component with a size 0(\B(N)\) whp (i.e., with 
probability tending to one as N —>■ oo), but for any p < p cr (c) the size of the largest connected 
component in Gjy(p,c) whp is at most 0(log \B(N)\). □ 

Hence, Theorem 11.11 may also tell us something about the "distances" between the com- 
ponents of a random graph when it is considered on the vertices of Z d . 

It is worth mentioning that the symmetry between c cr and p cr is most spectacular in the 
dimension one case, when p c (X) = 1- Notice, that d = 1 case is exactly solvable, and this is 
the only case when we know the distribution of |C|: 

P{|C| = k} = (1 -p) 2 kp k ~\ k>l. (1.7) 

Hence, if d = 1 we compute for all < p < 1 = p c (l) 

c cr {p) = (1.8) 



l+p 
1-c 



which also yields 

P cr (c) , 

1 + c 

for all < c < 1. (For more details on d — 1 case we refer to [ 1 1 J . ) 

Remark 1.1. For any fixed c function /3(p,c) is continuous atp = 0: if p = 0, i.e., when our 
graph is merely a classical G nfi / n random graph, then \C\ = 1 and ( ti.fij) becomes a well-known 
relation. 
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Furthermore, for any fixed c, if p = 0, it is not difficult to derive from fll.3p that 

a(0, c) = — *— r . (1.9) 

c - 1 + | logc| 

But log n/(c —1 + | logc|) is known (see Theorem 7a in [3]) to be the principal term in the 
asymptotics (in probability) of the largest connected component of G njC / n . This inevitably 
leads to the following (open) question. Will 

p{Ci(GV(p,c)) <« 1 log| J B(iV)|}^0 (1.10) 

hold also for all ol\ < a(p, c) when < p < p c and c < c cr (p)7 

It is easy to check that if c < c cr then the equation (II. 6p does not have a strictly positive 
solution, while (3 = is always a solution to (11.61) . This allows us to derive 

^ ' ciccr= 2 E^fy > °' (L11) 

which confirms that the phase transition remains to be of the second order for any p < p c , 
as it is for p = 0, i.e., in the case of classical random graph. 

More similarities and differences between our model and the "mean-field" case one can 
see in the following example. Let d — 1, in which case (11.7P holds. Introducing a random 
variable X with the first-success distribution 

P{X = k} = (l-p)p k - 1 , k = l,2,... , 

one can rewrite (jl.fip as follows (see the details in [TTj ) 

This equation looks somewhat similar to the equation obtained in [2] for the "volume" (the 
sum of degrees of the involved vertices) of the giant component in the graph with a given 
sequence of the expected degrees. Note, however, that in our model the critical mean degree 
when c = c cr and iV — >• oo is given according to (11.81) by 

2p + c cr = 2p + = 1 + (1.12) 

1+p 1+p 

which is strictly greater than 1 for all positive p < 1. This is in a contrast with the model 
studied in [2] , where the critical expected average degree is still 1 as in the classical random 
graph. 
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Although our model (it can be considered on a torus, in the limit the result is the same) is 
a perfectly homogeneous random graph, in the sense that the degree distribution is the same 
for any vertex, we study it via inhomogeneous random graphs, making use of the recently 
developed theory from [1]. The idea is the following. First, we consider the subgraph 
induced by the short-range edges, i.e., the edges which connect two neighbouring nodes with 
probability p. It is composed of the connected clusters (which may consist just of one single 
vertex) in B(N). Call a macro-vertex each of the connected components of this subgraph. 
We say that a macro-vertex is of type k, if k is the number of vertices in it. Conditionally 
on the set of macro-vertices, we consider a graph on these macro-vertices induced by the 
long-range connections. Two macro-vertices are said to be connected if there is at least one 
(long-range type) edge between two vertices belonging to different macro-vertices. Thus the 
probability of an edge between two macro- vertices and Vj of types x and y correspondingly, 



Below we argue that this model fits the conditions of a general inhomogeneous graph model 
defined in [1], find the critical parameters and characteristics for the graph on macro- vertices, 
and then we turn back to the original model. We use essentially the results from [1] to derive 
(11. 6p . The result in the subcritical phase (part i) of Theorem 11.11) does not follow by the 
theory in [I]; we discuss this in the end of Section 2.4. 

Notice also that in order to analyze the introduced model, we derive here some result on 
the joint distribution of the sizes of clusters in the percolation model (see Lemma I2TT1 below) . 
which may be of interest on its own. 

The principle of treating some local structures in a graph as new vertices ("macro- 
vertices"), and then considering a graph induced by the original model on these vertices 
appears to be rather general. For example, in [5] a different graph model was also put into 
a framework of inhomogeneous graphs theory by certain restructuring. This method should 
be useful for analysis of a broad class of complex structures, whenever one can identify local 
and global connections. Some examples of such models one can find in [7j. 

Finally we comment that our result should help to study a model for the propagation 
of the neuronal activity introduced in [12]. Here we show that a giant component in the 
graph can emerge from two sources, none of which can be neglected, but each of which may 
be in the subcritical phase, i.e., even when both p < p c and c < 1. In particular, for any 
< c < 1 we can find p < p c which allows with a positive probability the propagation of 
impulses through the large part of the network due to the local activity. 



is 
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2 Proof 



2.1 Random graph on macro- vert ices. 

Consider now the subgraph on B(N) induced by the short-range edges only, which is a purely 
bond percolation model. By the construction this subgraph, call it G^{p), is composed of a 
random number of clusters (of connected vertices) of random sizes. We call here the size of 
a cluster the number of its vertices (it may be just one). We recall here more results from 
percolation theory which we shall use later on. 

Let K N denote the number of the connected components (clusters) in (p), and let 

X = {X 1 ,X 2 ,...,X Kn } (2.1) 

denote the collection of all connected clusters Xj in G^(p). We shall also use Xj to denote 
the set of vertices in the i— th cluster. By this definition Ylf=i 1^1 = 1-^(^)1- 
Theorem [[1], (4.2) Theorem, p. 77] 

wh - K(p) : = (2 ' 2) 

a.s. and in L l as N —>■ oo. □ 

Note (see, e.g., [I]) that n{p) is strictly positive and finite for all < p < p c . Furthermore, 
in [T4] the large deviations property of is formulated as follows 



Theorem [[14J, Theorem 2] Given n(p) > e > 0, there exist <jj(e,p) > for j = 1,2 such 
that 

lim log P ( ^ > k(p) + e J = ai(e, p) 

and 



B n) \ B n 



lim log P ( ^ < K(p) — e ) = cr 2 (e, p) . 



n^oo 



B(n)\ \\B(n 



This theorem implies that for all < 5 < n(p) and all large iV there is a positive constant 
a = o-(5,p) such that 

K N 



n(p) 



\B(N)\ 

Define for any k > 1 and x > an indicator function: 



>5 \< e -\B(N)\_ (2.3) 



h.(x) 



1, if a; = fc, 
0, otherwise. 
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Proposition 2.1. For any fixed k > 1 



sl^'^^^' (2 - 4) 



o.s. and m L 1 as N oo. 



Proof. Let C(z), 2 G B(N), denote a connected (in G^\p)) cluster which contains vertex 
z. Then we write 



By the ergodic theorem 



K ^ *vi »i; - ^ k \ B ( N \\ 

JV i=l 2eB(iV) 



T^KT £ 4(|^)|) -P{|C| = fc} (2.6) 
a.s. as iV — ► 00. This in turn implies that convergence (I2.6P holds in L 1 as well, since 

°^jwm E 4(icwi) < 1. 

Hence, statement fl23| follows by ([23]) and Q. □ 
Given a collection of clusters X defined in (12. ip . we introduce another graph Gtv(X, p, c) 
as follows. The set of vertices of Gn(X.,p, c) we denote {t>i, . . . , f^}. Each vertex Vi is said 
to be of type which means that t>j corresponds to the set of \X^\ connected vertices in 
B(N). We shall also call any vertex V{ of G^(K,p, c) a macro-vertex, and write sometimes 

Vi = X { . (2.7) 

With this notation the type of a macro- vertex V{ is simply the cardinality of set Vi = JQ. The 
space of the types of macro-vertices is S = {1,2, . . .}. According to (12 ,4p the distribution 
of type of a (macro-) vertex in graph Gn(K,p, c) converges to measure // on S. The edges 
between the vertices of Gn(X.,p, c) are presented independently with probabilities induced 
by the original graph Gn(p,c). More precisely, the probability of an edge between any two 
vertices and Vj of types x and y correspondingly, is p xy (N) introduced in (11.13j) . Clearly, 
this construction provides a one-to-one correspondence between the connected components 
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in the graphs Gn(X.,p, c) and Gn(p,c): the number of the connected components is the 
same for both graphs, as well as the number of the involved vertices from B(N) in two 
corresponding components. In fact, considering conditionally on X graph Gjv(X,p, c) we 
neglect only those long-range edges from Gn(p,c), which connect vertices within each Vi, 
i.e., the vertices which are already connected through the short-range edges. 
Consider now 

*™ = l -( 1 -\imY =:! lSM- (2 ' 8) 

Observe that if x(N) — > x and y(N) — > y then 

K' N (x(N),y(N))^cxy (2.9) 

for all x, y G S. In order to place our model into the framework of the inhomogeneous 
random graphs from [1] let us introduce another (random) kernel 

K KN (x,y) = —2L- K? N [x,y), 

so that we can rewrite the probability p xy (N) in a graph Gn(K,p, c) taking into account the 
size of the graph: 

Pxy{N) = . (2.10) 

K N 

(We use notations from [1] whenever it is appropriate.) According to (12.21) and (12. 9ft . if 
x(N) — > x and y(N) — ► y then 

K KN (x(N),y(N)) ^ K (x,y) := CK(p)xy (2.11) 

as — > oo for all x, y G S . 

Hence, in view of Proposition 12.11 we conclude that conditionally on = t(N), where 
t(N)/\B(N)\ — > E(|C| _1 ), our model falls into the so-called "rank 1 case" of the general 
inhomogeneous random graph model G (t(N), Kt(N)) with a vertex space 

V = (S, fl, (v h . . . , V t (N))N>l) 

(see [I], Chapter 16.4). Note, that according to ( 11. II) function ji{k) (defined in ( 12.41) ) decays 
exponentially, which implies 

k G L l (S x S,fi x /i). (2.12) 
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Furthermore, it is not difficult to verify with a help of (12.21) and Proposition 12.11 that for any 
t(N) such that t(N)/\B(N)\ -> EdC]- 1 ) 



j j oo oo 

— -E{e(G^(X, p, c))\K N = t(N)} - - ^ ^ K ( x > vM^Mv), (2-13) 

^ ' j/=l cc=l 

where e(G) denotes the number of edges in a graph G. According to Definition 2.7 from pQ, 
under the conditions (I2.13p . (12.121) and (12.111) the sequence of kernels «t(jv) (on the countable 
space S x S 1 ) is called graphical on V with limit k. 

2.2 A branching process related to Gn(X-,p,c). 

Here we closely follow the approach from [1] . We shall use a well-known technique of branch- 
ing processes to reveal the connected component in graph Gn(X,p, c). Recall first the usual 
algorithm of finding a connected component. Conditionally on the set of macro- vertices, 
take any vertex u f to be the root. Find all the vertices {v\,v} 2 , ■ ■■,vj} connected to this 
vertex Vi in the graph Gn(K,p, c), call them the first generation of Vi, and then mark Vi as 
"saturated". Then for each non-saturated but already revealed vertex, we find all the ver- 
tices connected to them but which have not been used previously. We continue this process 
until we end up with a tree of saturated vertices. 

Denote tn(x) the set of the macro- vertices in the tree constructed according to the above 
algorithm with the root at a vertex of type x. 

It is plausible to think (and in our case it is correct, as will be seen below) that this 
algorithm with a high probability as iV — > oo reveals a tree of the offspring of the following 
multi-type Galton- Watson process with type space S = {1,2, . . .}: at any step, a particle 
of type x G S is replaced in the next generation by a set of particles where the number of 
particles of type y has a Poisson distribution Po(k(x, y)jj,{y)). Let p{x) denote the probability 
that a particle of type x produces an infinite population. 

Proposition 2.2. The function p(x), x G S , is the maximum solution to 

p{x) = 1 - e" E "=! >*>>vMv)pb) . (2.14) 

Proof. We have 

OO j 

V k(x, y)p(y) = cE(|C| _1 )x— — — - = cx < oo for any x, 

which together with ( 12.121) verifies that the conditions of Theorem 6.1 from [1] are satisfied, 
and the result ( 12. 14ft follows by this theorem. □ 
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Notice that it also follows by the same Theorem 6.1 from [T] that p(x) > for all x G S 
if and only if 

cB(\Cn f>Vy) = cE(\Cn f> 2 — ±— P{|C<l=y} = c E|C| > 1; (2.15) 

j/=i j/=i y\ \ ) y 

otherwise, p(x) = for all x G S. Hence, formula (11.21) for the critical value c cr (p) follows 
from fl2ZE5p . 

As we showed above, conditionally on so that Kn/\B(N)\ — > E(|C| _1 ), the sequence 
kk n is graphical on V. Hence, the conditions of Theorem 3.1 from [1] are satisfied and we 
derive (first, conditionally on K^, and therefore unconditionally) that 

C^GnQL^c)) p 
K N 

where p = YlT=i P( X )M X )- This together with (12.21) implies 

\B{Nj\ >E{l ° l )p - (216) 

Notice that here Ci(Gn(X.,p,c)) is the number of macro-vertices in the largest connected 
component of Gn(X-,P, c). 

2.3 On the distribution of types of vertices in Gn(K,p,c). 
Given a collection of clusters X (see (I2.1jl ) we define for all 1 < k < \B(N)\ 

K N 

A4=A4(X) = ^4(|X,|). 
1=1 

In words, A4 is the number of (macro-)vertices of type k in the set of vertices of graph 
Gjv(X,p, c). We shall prove here a useful result on the distribution of Af = (A/i, . . . ,Afx N )- 

Lemma 2.1. Set 

oo 

m = Y,*{\c\=n} 

n=k 

and fix v > 2 arbitrarily. Then for any fixed e > 



pjlM/ifiv-M^I > ek v Ji(k) for some 1 < fc < |S(JV)|} = o(l) (2.17) 



as iV — > oo. 
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Proof. Let us fix e > arbitrarily. Define a constant L so that eLq = E([C| l ). Then for 
all k > L 

ek u Jl(k)>i2(k), (2.18) 

and for any L > L Q 

P{\M k /K N - fi(k)\> ek v Jl(k) for some 1 < k < \B(N)\} (2.19) 
< P{ max \XA > L] 

l<i<K N 

+P{\M k /K N - /i(jfe)| > ek u Ji(k) for some 1 < k < L}. 

We shall choose later on an appropriate L = L(N) so that we will be able to bound from 
above by o(l) (as N — > oo) each of the summands on the right in (12.191) . 
First we derive 

P{ max |X,| > L} = P{ max |<7(*)| > L} < |S(JV)| P{|C| > L}, (2.20) 

1<i<Xjv zeB(N) 

where C(z) is an open cluster containing z. For a further reference we note here, that 
according to (II. ip for any < a < ((p) there is constant b > such that 

P{\C\ > L} < be- aL (2.21) 

for all L > 1, which together with (12.201) implies, in particular, that 

as N — > oo. 

Now we consider the last term in (I2.19p . Let us define for any < 5 < E(|C| _1 ) an event 



A 



K N 



Edcr 1 ; 



< 8 \ . (2.23) 



\B{N)\ 

Recall that according to (12.31) 

P(Aiv) > 1 - e-^ Bm = 1 - o(l) (2.24) 
as AT — > oo. Then we can bound the last term in (12.191) as follows 

p{\Af k /K N - ft(k)\ > ek u Ji(k) for some 1 < k < l} 

< p{\N k /K N - n(k)\ > ek u Jl(k) for some 1 < k < L } 

+ P| (U k /K N > e k v Ji(k) + fi(k) for some L < k < L j n A,iv} + P{A S , N } 

< p| (j\f k /K N > ek v Ji(k) + /i(k) for some L < k < L\ n A,at} + o(l), (2.25) 
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as N — > oo where the last inequality follows by Proposition 12.11 and bound (12.241) . Write 

P(k) :=p| (J^-fi(k)>ekrji(k) \ nA s , N y (2.26) 
Clearly, we have by (12. 25ft : 

L 

pUAf k /K N -fi(k)\ > ek v Jl(k) for some 1 < k < l| < ^ P ( k ) + o(l), (2.27) 

k=L +l 

as N — > oo. Substituting now fl2~27j) and f[2~20l) into fl2~T9l) we derive 

P{\Af k /K N - fi(k)\ > ek v Ji(k) for some 1 < k < \B(N)\} (2.28) 

L L 

<|i?(iV)|P{|C|>L}+ P(k) + o(l)<\B(N)\m + E P (k) + o(l) 

k=L +l k=L +l 

as N — > oo. 

Next we shall find an upper bound for P{k). Due to the definition (12.231) of As ; n, we 
have 

P(k) < P |A4 > (k(p) - < y)|S(JV)|(eA; ,/ /i(A:) +/x(A;)) J. (2.29) 

We shall use the following special case of the Talagrand's inequality. 
Proposition 2.3. For every r£K and t > 

P{^<r- ( }p{^>r}<exp{-^}. (2.30) 

Proof. We shall derive this result as a corollary to the Talagrand's inequality [8] which we 
cite here from the book [6], p. 40. 

Theorem. [Talagrand's Inequality] Suppose that Zi,...,Z n are independent random vari- 
ables taking their values in some sets A\,...,A n , respectively. Suppose further that W = 
f(Zi,...,Z n ), where f : Ai x ... x A ra — > M. is a function such that, for some constants Ck, 
k = 1, ...,n, and some function ty, the following two conditions hold: 

1) If z, z' G A = Y\i Aj differ only in the i — th coordinate, then \ f(z) — f{z')\ < ci. 

2) If z G A and r G K with f(z) > r, then there exists a set J C {1, ...,n} with 
J2ie.J c "i — ^( r )' suc h that for ally G A with yi = z { when i G J, we have f(y) > r. 
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Then, for every reR and t > 0, 

P(W < r - t)P(W > r) < e^ 4 ^. (2.31) 

We shall show now that function M k satisfies the conditions of this theorem. Let 
{ei, . . . , e n } be the set of all edges from the lattice Z d which have both end points in B(N). 
Define 

_ _ f 1, if e { is open in G%(p), 
\ 0, if e, is closed in G s N (p) . 

According to the definition of our model, Zj G Be(p), % — 1, . . . , n, are independent random 
variables, and 

M k =Af k (Z u ...,Z n ) 

since A4 is the number of the components of size k (open fc-clusters) in G s N (p), which is 
defined completely by Z 1 , . . . , Z n . Furthermore, it is clear that removing or adding just one 
edge in G s N (p), may increase or decrease by at most one the number of fc-clusters. Hence, 
the first condition of the Talagrand's inequality is satisfied with c« = 1 for all 1 < i < n: if 
configurations z, z' G {0, l} n differ only in the i th coordinate, then 

\N k {z)-N k {z')\<l. 

Next we check that the second condition is fulfilled as well, and we shall determine the 
function \E'. Assume, z G {0, 1}™ corresponds such configuration of the edges in B(N) that 
Nk(z) > r , for some r > 1, i.e., there are at least r clusters of size k. Let {ej,j G J} C 
{ei, . . . , e n } be a set of edges which have at least one common vertex with a set of exactly 
r (arbitrarily chosen out of Afk(z)) clusters of size k. Clearly, |J| < 2dkr, and for any 
z' G {0, 1}™ with z'j = Zj if j G J, we have 

M k {z') > r, 

proving that the second condition of the Talagrand's inequality is satisfied as well with 
\l/(r) = 2dkr, since 

= \J\<2dkr. (2.32) 

Hence, the inequality fl2T30|) follows by fr2~3B . □ 

Set now 

k N = (K(p)-5)\B(N)\, 
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(2.33) 



and consider the inequality (I2.3(jp with 

r = k N (ek u Ji(k) + fJ,(k)\ 

r-t = k N (^k u Jl(k) + n(k)j. 
First we derive for any k > L (in which case ek u Jl(k) > fi(k)) 

p{*1 < r - t } > p{*1 < §*„„(*)} > 1 - pH). (2.34) 
where we used the Chebyshev's inequality. Recall that by Proposition 12.11 (and (12.21) ) 

Ml K ^ k) 

as iV — ► oo. Hence, choosing < 5 < n(p)/10 we have 

EA4 EA4 



|(k(p)-5)|S(JV)|/x(A;) 
for all large iV, which together with (12. 34ft implies 



< 3/4 



p|A4<r-t|>^ (2.35) 

for all large N. Using the last bound in the Talagrand's inequality (I2.30P with r and t defined 
in (I2.33p . we derive for all large N when k > Lq (and therefore ek u Jl(k) > fi(k)) 

Pk>r}<(pk<r-t}) -p{-dU (2 ' 36) 



dk N k»]i(k)y 



dk [k N [ek u + fi(k) 



< 4 exp 



, 4CTp {_ £ ^} = 4exp {_ £ (^) |B(JV)| ^, w 

Substituting (12T36D into (12T29D we get 

P(k) < 4exp {-a|£(JV)|fc"- 1 Jl(k)} , 

14 



where 

e{k{p) - 5) 



a :-- 



64d 

The last bound combined with (12.281) yields 



p{\Af k /K N -n(k)\ > ek v Jl(k) for some 1 < k < \B(N)\ } (2.37) 

L 

< \B(N)\Ji(L)+AJ2 exp {-alBiN)^ 1 ]l(k)} +o(l), 

k=L 

as N — > oo for any L > L . 

Next we shall show that for any A > one can choose a finite constant L Q and numbers 
L = L(N) such that 

L(7V) 

^ exp{-a|S(A^)|r- 1 /I(A;)} < A, (2.38) 

k=L 

for all large N, and 

|S(iV)|/I(L(iV))^0, (2.39) 

as N — ► oo. This together with (12.371) will clearly imply the statement of Lemma. 
We claim that both (TOHj) and (|2T39|) hold with 

L(N) = min j* : < j^^} , (2-40) 

where 

v-2 
a = — 

is positive by the assumption of Lemma. Observe, that k a Jl(k) — > 0, as k — > oo for any fixed 
a due to the exponential decay (II. ip . This yields that L(N) — > oo as iV — > oo, which in turn 
implies that there exists 

lim J1(L(N))\B(N)\ < lim L(N)~ a = 0, 

AT— »oo N^oo 

and (1239]) follows. 

To prove (T2~38D first we note that by the definition (I2T401 of L(JV) 

- l) a /I(L(iV) - 1) > —Ljj. (2.41) 



15 



Recall that according to Lemma 6.102 from [4] (p. 139), for all n, m > 

' P(|C| =n + m) >p(l -pY 2 — P(|C| = m)-P(\C\ =n). (2.42) 



m + n m n 

When m = 1 the inequality (12.421) implies 

P(|C| = n + 1) > p(l - p) 2(d - 1) P(|C| = n), (2.43) 

for all n > 0. This clearly yields 

£(L(JV)) > p{\ - p) 2 ^Jl(L(N) - 1). (2.44) 

Notice that 7 := p{\ - p) 2 ^-^ < p < 1 for all d > 1. Combining (EH]) with (EM) we 



immediately get 



L(N) a Ji(L(N)) > j^-, (2.45) 



and also by the definition (12.401) for all k < L(N) 



Making use of (I2.45P and (I2.46[) we derive 

L(N) 

exp{-a|5(iV)|A; I/ - 1 /I(A;)} (2.47) 

k=Lo 

L(N) 

< ex P {-ajk"- 1 - } < ai exp {-a 7 L^ Q } , 

k=Lo 

where a\ is some positive constant independent of Lq. It is clear now, that for any A > 
we can fix L so large that (I2.47P would imply (I2.38p . and in the same time L will satisfy 
(I2.18P and L < L(N). This completes the proof of the lemma. □ 

2.4 Proof of Theorem 11.11 in the subcritical case c < c cr (p) . 

Let us fix < p < p c and then c < c cr (p) arbitrarily. Given X let again V{ denote the macro- 
vertices with types i — 1,2, ... , respectively, and let L denote a connected component 
in Gn(X,P, c). Consider now for any positive constant a and a function w = w(N) 

P jCi (G N (pc) j > au>j = P < max Y Xi > aw \ . (2.48) 
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We know already from (I2.16P that in the sub-critical case the size (the number of macro- 
vertices) of any L is o(N) with probability tending to one as iV — > oo. Note that when 
the kernel k(x, y) is not bounded uniformly in both arguments, which is our case, it is not 
granted that the largest component in the subcritical case is at most of order log|JB(JV)| 
(see, e.g., discussion of Theorem 3.1 in [1]). Therefore first we shall prove the following 
intermediate result. 

Lemma 2.2. If c < c cr (p then 

p{c 1 (G iV (X,p,c)) > ^(iV)! 1 / 2 } =o(l), (2.49) 

as N — > oo. 

Proof. Let us fix e > and 5 > arbitrarily and introduce the following event 

B N = A s , N n( max \X t \ < -^log \B(N)\ ] (2.50) 

\l<i<K N ((p) J 



< ek v jl(k) 



K N 

According to fl2~24l) . fT2T22|) and (ETTTl) we have 

P{B N } = l-o(l) (2.51) 

as N — > oo. 

Recall that tat(o;) denote the set of the macro- vertices in the tree constructed according 
to the algorithm of revealing of connected component described above. Let |tjv(x)| denote 
the number of macro- vert ices in tn(x). Then we easily derive 

p{c 1 (g w (X,j9,c)) > Ifi(^)! 1 / 2 } <p| i max^|r JV (|X i |)| > ^(iV)) 1 / 2 | B N \ + o(l) 

(2.52) 

\B(N)\ 

< \B(N)\ (6 + E(|C|- 1 )) (Kk) + ek»il{k))P{\T N (k)\ > ^(iV)) 1 / 2 | B N } + o(l) 

fc=i 

as N — >• oo. We shall use the multi-type branching process introduced above (Section 2.2) 
to approximate the distribution of |rjv(fc)|. Let further X c,p (k) denote the total number of 
the particles (including the initial one) produced by the branching process starting with a 
single particle of type k. Observe that at each step of the exploration algorithm, the number 
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of new neighbours of x of type y has a binomial distribution Bin(N y ,p xy (N)) where N' is 
the number of remaining vertices of type y, so that N' y < M y . 

We shall explore the following obvious relation between the Poisson and the binomial 
distributions. Let Y UiP e Bin(n,p) and Z a G Po(a), where < p < 1/4 and a > 0. Then for 
all k > 

P{Y n , p = k}<(l + Cp 2 ) n P{Z nj z_ = k}, (2.53) 

where C is some positive constant (independent of n, k and p). Notice that for all 1 < x, y < 
C§)log|5(iV)| 

(\ x v 
l ~wm =WN)\ xy{l+0{l)) ^ (2 - 54) 

and clearly, p xy (N) < 1/4 for all large N. Therefore for any fixed positive E\ we can choose 
small e and 5 in (12.501) so that conditionally on Bn we have 

N' P _ xy ~ N L < (My) + y v eim)<x, v) (2.55) 

1 Pxy{^ ) 

for all large N. 

We shall use the following property of measure 

u(k) - 1 P ^l = fc > 
MW - E(|C|-i) k 

defined in Proposition 12.11 Recall, that along with the result fll.il) it is also proved in [4] 
that for all < p < p c 

C(p) = i«m fW00 f-^logP{|C| >n}^j . (2.56) 

Hence, (11.11) and (12.561) immediately imply the existence and equality of the following limits 
for all < p < p c 

C(p) = Urrin^oo f-ilog//(n)^ = lim^^ (~ log/2(ra)^ , (2.57) 

i.e., that both /i(n) and /I(n) decay exponentially fast, and moreover with the same exponent 
in the limit. Let us write further 

Ms/) = Vp(y)> My) = h(y), y) = y) = cn(p)xy, 
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emphasizing dependence on p and c. The result ( I2.57P allows us to choose for any positive 
e 2 and p < p' < p c a positive S\ = Ei(e 2 ,p') such that 

H P (y) + y v e x Ji p {y) < (1 + e 2 )^{y). (2.58) 

Setting now d := (l+e 2 )^k c we derive from (12.551) with a help of (I2.58p . that conditionally 
on Bjy with an appropriate choice of constants 

Px (N) 

N ' y l-p (N) ~ ^ 1 + £2 ^p'^^p( x ^) = Vp'(y)Kc>,p>{x,y)- (2.59) 

Recall that above we fixed p and c < c cr (p), where c cr (p) is strictly decreasing and continuous 
in p. Furthermore, function n(p) is analytic on [0,p c ). Hence, we can choose p' > p and 
d := (1 + e 2 ) so that 

c < d < c cr {p') < c cr (p), (2.60) 

and moreover d and p' can be chosen arbitrarily close to c and p, respectively. Now according 
to fl233|) and (12391) 

p { Y N>,p xy (N) >k}<(l + Cp xy (N) 2 fy P{Z ~ Pxy(N) > k} (2.61) 

V l-Pxy(N) 

< (l+Cp xy (N) 2 )\ B ^P{Z, pl{y)Kc/y(x ^ > k}. 

Hence, if conditionally on Bn at each (of at most \B(N)\) step of the exploration algorithm 
which reveals T^(k), we replace the Bin(N y ,p xy (N)) variable with the Po (/i p i(y)K c ' tP '(x,y)) 
one, we arrive at the following bound using branching process: 

P{\r N (k)\>\B(N)\ 1 / 2 \B N } (2-62) 

2 \ \B(N)\ 2 



< 1 + C 



1/2 



This together with (I2.54p implies 

P {Mife) I > ^(iV)! 1 / 2 I B N ) < e ^°9\B{N)\) A P ^x c '*\k) > ^(iV)! 1 / 2 } , (2.63) 

where b\ is some positive constant. Substituting the last bound into (12.521) we derive with a 
help of (12381) 

p{c 1 (G A r(X,p,c)) > I^KiV)! 1 / 2 } (2.64) 
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\B(N)\ 

< b 2 \B(N)\e b ^ B ( N ^ kn P >(k)P{x c '' p '(k) > |5(A0r /2 j + o(l) 

fc=i 

as N —>■ oo, where b 2 is some positive constant. By the Markov's inequality 

P{X c '*'(k) > ^(iV)! 1 / 2 } < 2 -l^)l 1/2 E ^ c ' y W (2.65) 
for all z > 1. Denote /i z (£;) = Ez* c ' P then with a help of (12.651) we get from (12.641) 

p{c 1 (g jv (X,p,c)) > ^(iV)! 1 / 2 } (2.66) 

|B(JV)| 

<b 2 \B{N)\e b ^\ B ^\~\ B ^\ 1/2 k^{k)h z {k) + o{l). 

k=l 

Now we will show that there exists z > 1 such that the series 

00 00 
D z (c',p') = Y k ^P'(k)h z (k) = KfrT^VAM = k}h z (k) (2.67) 

k=l k=l 

converge. This together with (I2.66P will clearly imply the statement of the lemma. 

Note that function h z (k) (as a generating function for a branching process) satisfies the 
following equation 

h z (k) = zexp {J27=i Kc',p>(k,x)n P >(x)(h z (x) -1)} 

= zexp{c f K(p')k(Y^^ =1 xii p '(x)h z (x) - k(p'Y 1 )} 

= zexp {c'K(p')k(B z (c',p') - nip')- 1 )} . 
Multiplying both sides by k(j, p >(k) and summing up over k we find 

D z (c',p') = Y k ^(k)zexp {c'K(p')k(D z (c',p') - ^p')" 1 )} 

k=l 

oo 

= J2 <PT l Vp>{\C\ = k}zexp {c'K(p')k(D z (c',p') - ^p')- 1 )} , 

k=l 

where we also used the definition of fi(k). Let us write for simplicity D z = D z (c,p). Hence, 
as long as D z is finite, it should satisfy equation 

D z = n(p)- l zVe c \ c \ {K{p)D *- l \ (2.68) 
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which implies in turn that D z is finite for some z > 1 if and only if (I2.68P has at least one 
solution (for the same value of z). Notice that by the definition (12.671) 

D Z >D 1 = ^p)- 1 = (EflCI- 1 ))- 1 (2.69) 

for z > 1. Let us fix z > 1 and consider equation 

y/z = Ee c|c|(y - 1} =: F(y) (2.70) 

for y > 1. Using the property (11. ip of the distribution of |C| it is easy to derive that function 
F(y) is defined on [0, ({p)/c) where it is finite, increasing and has positive second derivative. 
Compute now 

§ j F(y)\ y=1 = cE\C\ = ±. (2.71) 

Hence, if c < c cr then there exists z > 1 such that there is a finite solution y to (12.701) . and 
therefore (12.681) also has at least one solution for some z > 1. Taking into account condition 
(I2.60p . we find that D z (c',p') is also finite for some z > 1, which finishes the proof of the 
lemma. □ 

Now we are ready to complete the proof of (11.41) . following almost the same arguments 
as in the proof of the previous lemma. Let Sn(x) = J2xi£r N (x) denote the number of 
vertices from B(N) which compose the macro- vert ices of tn(x). Denote 

B' N := B N n (d (d N (X,p,cf) < |£?(A0| 1/2 ) • 
According to (12.511) and Lemma 12.21 we have 

P{B' N } = l-o(l). 

This allows us to derive from (I2.48P 

P {Ci(Gtf(p,c)) > aw} < P ^max S N (Xi) > aw \ B' N j + o(l) (2.72) 

\B(N)\ 

< \B(N)\(6 + EdC)- 1 )) + ek v Ji(k))¥ {S N {k) > aw | B' N } + o(l). 

k=l 

Let now S c,p (y) denote the sum of types (including the one of the initial particle) in the 
total progeny of the introduced above branching process starting with initial particle of type 
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y. Repeating the same argument which led to (I2.62p . we get the following bound using the 
introduced branching process: 

P {S N (k) > aw | B' N ] 

/ 2 s b!\B(N)\^\B(N)\ 

< 1 1 + C ( max p xv (N)) ) P \s c ' ;p '(k) > aw] 

\ \x,y<2\og\B(N)\/ap) J > 



as N — > oo, where we take into account that we can perform at most y^\B(N)\ steps of 
exploration (the maximal possible number of macro- vertices in any L conditioned on B' N ). 
This together with (I2.54p implies 

P {r N (k) > aw | B' N } < (1 + o(l))P {s c '' p '{k) > aw} (2.73) 

as N — > oo. Substituting the last bound into (I2.72p we derive 

\B(N)\ 

PWi\G N (p,cj) > aw} < b\B(N)\ kfi p ,(k)P {S c '' p \k) > aw} + o(l) (2.74) 

fe=i 

as iV oo, where 6 is some positive constant. Denote g z (k) = ~Ez sc ' P then similar to 
( 12361) we derive from (J2Z3D 

|B(JV)I 

pjc^Gjvfec)) >aw(iV)} < 6|B(iV)| ^ k[i p i(k)g z (k)z~ aw{ - N ^ + o(l). (2.75) 

fe=i 

We shall search for all z > 1 for which the series 

oo 1 oo 

A z (c',p') = £>M*)&(*) = -JlsJ2 P P'i\ C \ = k ^( k ) 
k=l > k=l 

converge. Function g z (k) (as a generating function for a certain branching process) satisfies 
the following equation 

g z (k) = z k exp {J27=l K c\ P '(k } x)/i p >(x)(g z (x) - 1)} 

= z k exp{c , K(p')kJ2™=i x l J 'p'( x )(gz( x ) ~ !)} 
= z k exp {c'k(K(p')A z (c' ,p') - 1)}. 
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Multiplying both sides by kn P >(k) and summing up over k we find 

oo 

A z (c',p') = J2^ P ik)z k exp{c'k(K(p')A z (c',p') - 1)}. 

k=l 

Denoting for simplicity A z = A z (c,p), we can rewrite the last equation as follows: 

oo 

A z = kfi P (k)z k exp {ck{K{p)A z - 1)} . (2.76) 
fe=i 

It follows from here (and the fact that A z > 1/k(p) = A\ for all z > 1) that if there exists 
z > 1 for which the series A z converge, it should satisfy by fll.il) 

z < e c{p) . (2.77) 

According to (12.761) . as long as A z is finite it satisfies the equation 

A z = (k( P ))- 1 E (z^ e «W(*(pM*-i)) , 

which implies that A z is finite for some z > 1 if and only if the last equation has at least 
one solution 

A z > A, = l/«(p). (2.78) 
Let us fix z > 1 and consider equation 

y = E(z^e c \ c ^) (2.79) 

for y > 1. It is easy to check that at least for some y > 1 and z > 1 function 

f(y,z) :=E(zWe^y-V) 

is increasing, it has all the derivatives of the second order, and ^f(y, z) > 0. Compute now 

^-f(y,z)\ y=1 , z=1 = cE\C\ = ^. (2.80) 
ay c cr 

Hence, if c > c cr there is no solution y > 1 to (I2.79P for any z > 1. On the other hand, if 
c < c cr then there exists 1 < Zq < e^> such that for all 1 < z < z there is a finite solution 
y > 1 to (12.791) . We shall find z as the (unique!) value for which function y is tangent to 
f(y,Zo) if 2/ > 1- 
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First we rewrite (12.791) as follows. Set 

a = - log z, 

c 

then (12.791) is equivalent to 

y = Ee c\C\(y-l+a)^ ^ 281 ^ 

which after the change x = y — 1 + a becomes 

x + l-a = Ee clclx . (2.82) 

Here on the right we have a convex function with a positive second derivative (for all x < 
C(p)/c). Notice also that by the assumption 

^-Ee c|c|x L =0 = EclCI < 1. 
ox 

Hence, there exists unique y > such that 

— Ee c ^ x \ x=yo = E (c\C\e c \ c ^ = 1. (2.83) 

Define now 

a :=l+y -Ee c ^ yo , (2.84) 

which is strictly positive due to the preceeding argument. Clearly, function x + 1 — is 
tangent to Ee^K Hence, for all a < ao equation ( 12. 82ft has at least one solution, which 
implies due to (12.811) that for all 

z < z := e mo = exp{c(l + y -Ee c|c|w )} (2.85) 

equation (I2.79P has also at least one finite solution y > 1. This yields in turn that A z is 
finite for all z < z . 

Now taking into account that d > c and p' > p can be chosen arbitrarily close to c and 
p, respectively, we derive from (12.751) that for all 1 < z < z 

p{Ci(Gjv(p,c)) >aw{N)} < b(z)\B(N)\z- aw{N) + o(l) (2.86) 

as iV — > oo, where b(z) < oo. This implies that for any 

a > 1/logzo = (c + q/ - E ce c|c|l ' )~ 1 
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and w(N) = log \B(N)\ 



p{Ci(G;v(p,c)) >alog|fl(iV)|} = o(l) (2.87) 
as A" — > oo, which proves (11.41) . □ 

To conclude this section we comment on the methods used here. It is shown in [9] that 
in the subcritical case of the classical random graph model G n ^ c / n (i.e., p = in terms of our 
model) the same method of generating functions leads to a constant which is exactly a(0, c) 
(see ( II. 9p ). The last constant is known to be the principal term for the asymptotics of the 
size of the largest component (scaled to log A/") in the subcritical case. This gives us hope 
that the constant a(p, c) is close to the optimal one also for p > 0. 

Similar methods were used in [10] for some class of inhomogeneous random graphs, and 
in pQ for a general class of models. Note, however, some difference with the approach in pQ. 
It is assumed in pQ, Section 12, that the generating function for the corresponding branching 
process with the initial state k (e.g., our function g z (k), k > 1) is bounded uniformly in k. 
As we prove here this condition is not always necessary: we need only convergence of the 
series A z , while g z {k) is unbounded in k in our case. Furthermore, our approach allows one 
to construct constant a(p, c) as a function of the parameters of the model. 

2.5 Proof of Theorem 11.11 in the supercritical case. 

Let Ck denote the set of vertices in the fc-th largest component in graph Gn(p,c), and 
conditionally on X let Ck denote the set of macro- vert ices in the fc-th largest component 
in graph Gat(X,p, c) (ordered in any way if there are ties). Let also Ck and Ck denote 
correspondingly, their sizes. According to our construction for any connected component L 
in Gn(X.,p, c) there is a unique component L in Gn(p,c) such that they are composed of 
the same vertices from B(N), i.e., in the notations (12.71) 

L = u x l6 Z U.ex, {z} =: V(L). 

Next we prove that with a high probability the largest components in both graphs consist 
of the same vertices. 

Lemma 2.3. For any < p < p c if c> c cr (p) then 

P{d = V(Cx)} = 1 - o(l) (2.88) 

as N — » oo. 
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Proof. In a view of the argument preceeding this lemma we have 

P{d ^ V(Cx)} = P{Ci = V(C k ) for some k > 2}. 

According to Theorem 12.6 from [lj, conditions of which are satisfied here, in the su- 
percritical case conditionally on Kn such that Kn/\B(N)\ — > E(|C| _1 ), we have whp 
C 2 = 0(\og(K N )), which by Proposition 12.21 implies C 2 = 0(log \B(N)\) whp. Also we 
know already from (I2.16P that in the supercritical case C\ = 0(\B(N)\) whp, and therefore 
C\ = 0(\B(N)\) whp. Hence, for some positive constants a and b 

P{d ^ V(Ci)} = P{C X = V(C k ) for some k > 2} (2.89) 

< pjUax \V(C k )\ >b\B(N)\ S J n (m&xCk < a\og\B(N)\^ + o(l). 
It follows from fl2"^2l that 



P{ max \X t \> y/\BjN)\} = o(l) 

l<l<Kjv 



as N —>■ oo. Now we derive 



P jhoMK|V(C fc )| >b\B(N)\J n (maxC k < a\og\B(N)\^j\ (2.90) 

<pj (max \V{C k )\ > b\B(N)\ J n (maxC fc < alog\B(N)\ ) D ( max |JQ| < (JV) | 

+o(l) < P {y/\B(N)\a\og |S(JV)| > 6|S(JV)|} + o(l) = o(l). 
Substituting this bound into (I2.89j) we immediately get (12.881) . □ 

Conditionally on C\ = V(C\) we have 

\B(N)\ = \B(N)\ EiJl |Aj|l{Aj eCi} 

= jnm EE? 01 = G Ci} (2.91) 

= ji^ k£#{Xi G & : |X<| = k}. 

Note that Theorem 9.10 from p] (together with (12.21) in our case) implies that 

Mk) ■= 4-#{ X * G : = fc > ^ P(%(*0 (2-92) 
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for each k > 1, where p(k) is the maximal solution to (12.141) . 
We shall prove below that also 

|B(AT)| oo 

W N := ku N (k)^^2kp(k)fi(k)=:l3 (Edq- 1 )) -1 . (2.93) 

k=l k=l 

Observe that according to (I2.14p constant (3 (defined above) is the maximal solution to 

oo oo 

/3 = E(|C|- 1 )^A;p(A;)/i(A;) = E(|C| _1 ) ^ k ^1 - e ~ W^Mj/)p(^ fi{k) 

k=l k=l 



This proves that (3 is the maximal root of ( 11. 61) . Then ( 12.93[) together with ( 12.2ft will allow 
us to derive from ( I2.9ip that for any positive e 

cJg k ( p ,c) 



p {| \B(N)\ ' -!>\><\ey = nty^ 

as N — > oo. This combined with Lemma [2.31 would immediately imply 

Ci(g n (p, c 



\B(N)\ 



P - P, (2-94) 



and hence the statement of the theorem follows. 

Now we are left with proving (12.931) . For any 1 < R < \B{N)\ write W N := W§ + w§, 
where 

R \B(N)\ 
W$ := ^ku N (k), w§:= k "N(k). 

k=l k=R+l 

By (12.921) we have for any fixed R > 1 

^iv * ( 2 - 95 ) 
fc=i 

as iV — ► oo. Consider w§. Note that for any k > 1 

E^(A;)<E^^/ fc (|X J |) = E^^i-l- £ J*(|C(*)|), (2.96) 



27 



S,N j 



where C(z) denotes again a connected cluster in the bond percolation model on B(N) with 
a probability p of bound. Using events As,n together with bound (12.241) . we obtain from 
fl^96|) for any fixed < 5 < E(|C|- 1 )/2 and k > 1 

< v(\C\\-6 P {I C I = *} + IBWIPfA^}. 
Bound (12.241) allows us to derive from here that 

Eu N (k) < A 1 (P{\C\ = k} + e - ailBm ) (2.97) 

for some positive constants A\ and a\ independent of k and N. This together with (11.11) 
yields 

\B(N)\ 

= kEu N (k) < A 2 e~ a2R (2.98) 

k=R+l 

for some positive constants A 2 and a 2 - 

Clearly, for any £>0we can choose i?o so that for all R > Rq 

oo 

kp(k)p(k) < e/3, 

fc=i?+i 

and then we have 

oo 

PII^tv-^M^M^I >4 (2-99) 



fc=i 

R 



= P{\(W*-Y,kp(k)p(k))+w*- kp(k)p(k)\>e} 

k=l k=R+l 
R 

< P{\W* - k P (k)p(k)\ > e/3} + P{w§ > e/3}. 

k=l 

Markov's inequality together with bound (12.981) gives us 

P K>E/3} <jM<^- a (2 . 100 ) 
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Making use of (12.1001) and (I2.95P we immediately derive from (12.991) 

P{\W N -Y^kp(k)n(k)\>e}<o(l) + — (2.101) 

fc=i £ 

as N — > oo. Hence, for any given positive e and Eq we can choose finite R so large that 

oo 

lim V{\W N - V kp(k)n(k)\ > e} < e . (2.102) 

N^oo ' J 

k=l 

This clearly proves statement (12.931) . and therefore finishes the proof of the theorem. □ 
Acknowledgment T.T. thanks MSRI for the hospitality at the beginning of this project. 
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